Optical response of a misaligned and suspended Fabry-Perot cavity 
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The response to a probe làser beam of a suspended, misaligned and detuned optical cavity is 
examined. A five degree of freedom model of the fluctuations of the longitudinal and transverse 
mirror coordinates is presented. Classical and quantum mechanical effects of radiation pressure are 
studied with the help of the optical stiffness coemcients and the signals provided by an FM sideband 
tcchnique and a quadrant detector, for genèric vàlues of the product vjt of the fluctuation frequency 
times the cavity round trip. A simplified version is presented for the case of small misalignments. 
Mechanical stability, mirror position entanglement and ponderomotive squeezing are accommodated 
in this model. Numerical plots refer to cavities under test at the so-called Pisa LF facility. 



I. INTRODUCTION AND NOTATION 

Optical cavities are generally studied by assuming a 
single mode excitation and ignoring the photon scatter- 
ing by mirror reflections into other modes. A single mode 
description is no more reliable when a slight misalignment 
is sufïicicnt to cxcitc difïcrcnt modes. This situation is 
met in almost concèntric or plane-parallel or confocal 
configurations. In other cases the weak amplitudes of 
modes falling outside the resonance bandwidth are of in- 
terest. For example, the sensitivity of intcrfcrometers 
with cavities placed in their arms depends on the contri- 
bution of higher modes as well as the error signals used 
for longitudinal and angular alignments. 

Optical cavities are generally stabilized, in length, by 
the Drever-Pound (D-P) technique Q working with odd 
harmonies of the phase modulated làser beam. 

In order to deal with a large variety of situations the 
model discussed in this paper accounts for misalignment, 
detuning and a genèric spectrum of harmonies. Faced 
with the possibility of working with approximate expres- 
sions it has been preferred to simplify exact solutions at 
the end of the calculation. This option avoids difficulty 
of making adequate approximations in presence of a large 
number of parameters. 

This strategy can be useful for the design studies nec- 
essary for the development of future gravitational anten- 
nas. It gives the possibility to investigate noise contri- 
butions coming from all optical and mechanical degrees 
of freedom. It can be also used for studying instabilitics, 
optical spring effect, entanglement and radiative pressure 
squeezing associated to both axial and angular fluctua- 
tions for any degree of detuning, misalignment and mis- 
match. 

The present work grew up from the study of short 
and large spot size resonators 0] implemented at LFF 
(Low Frequency Facility 0) a facility dedicated to test- 



ing new mechanical suspensions, controls and mirrors for 
the VIRGO interferometric gravitational antenna, and 
studying the effects of radiation pressure, mirror and sus- 
pensions thermal noise. 

Main features of the LFF are the use of suspended 
mirrors and the possibility of confining large section cav- 
ity modes. The mirrors hang from multipendula which 
guarantee a dràstic reduction of the seismic noise above 
the resonance freqüències of the mechanical modes. The 
phase-modulated light reflected by the cavity is used by a 
Pound-Drever apparatus Q both for stabilizing the cav- 
ity length, and measuring the noise spectrum. Several 
papers analyzed the dynamic and the alignment of cav- 
ities sharing some of the LFF features 0, TU |(| . Numer- 
ous specialized studies have been produced by research 
groups of VIRGO, LIGO, TAMA and GEO projeets Q. 

Suspended cavities have been analyzed by several au- 
thors in different contexts, all sharing the common fea- 
ture of using a system of Langevin equations for both the 
mechanical and electromagnètic modes. The coupling of 
cavity and mechanical modes is represented by suitable 
potentials 0, E} leading to a complex interplay between 
cavity mode amplitudes, mirror positions and orientation 
fluctuations. In this paper the resonator is regarded as 
a mechanical Langevin system driven by thermal sources 
and shot noise. This is done by-passing the Hamiltonian 
approach and hiding the optical modes fluctuations into 
the mechanical ones by generalizing an approach intro- 
duced in [ldj . So doing, the Langevin system contains 
ponderomotive terms, connected with the classical part 
of the làser beam and the shot noise. The seismic noise af- 
fecting the mechanical system has been neglected. Once 
known, its local spectral density can be easily added to 
the thermal one. 

Thermal motion of mechanical oscillators has been 
modelled as Standard Brownian motion |ïïj , possibly cor- 
rected by Diosi for preserving the quantum mechanical 
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FIG. 1: Typical optical layout of the apparatus examined 
in the present paper. The làser beam is phase modulated at 
frequency A by the modulator M (A) . The modulated bema 
drives the cavity. The front M\ and back M2 mirrors are 
suspended to multipendulum chains. The beam reflected by 
the front mirror is sent to a photodetector PD, and the pho- 
todetection signal is demodulated before going to a spectrum 
analyzer. The same signal is sent to a control system (not 
shown) which provides a feedback signal applied to mirror 
Mi. The noise of the system is studied by spectral analysis. 



commutation relations [Tal, or by non-Lindblad master 
equations (ME) 0, Il4j. Accordingly, in the present 
model different thermal correlation functions have been 
introduced. 

The quantum field fluctuations (shot-noise) are ac- 
counted for by splitting each mode amplitude in a classi- 
cal and a quantum parts 0, IÏ5| and relying on the input- 
output theory 16]. 

Radiation pressure can lead to mechanical instabilities, 
as predicted by Braginsky and Manukin 01 • Acting on 
the suspended mirrors it provides a spring action which 
either depresses or reinforces any perturbation |Tot Esl 
lï^ . l2(ï| . It may also be used to mechanically entangle 
the two mirrors [2Ï| or to enhance the squeezing of the 
output field [12. 

In Fig^ it is represented the typical optical layout of 
the apparatus that will be examined. Main features of 
the present model are: i) the multimode description of 
the cavity field; ii) the inclusion of radiation pressure and 
shot noise terms; iii) the description of suspensions and 
mirrors in terms of mechanical modes. 

Here and in the following J = 1,2 labels the mirror, 
x the axial and q — y, z the transverse coordinates. The 
analysis is focused on the fluctuations of the J's mirror 
orientation {59 j q ) and displacements (óxj, Sej q ), com- 
bined in the parameters 



Sípj = {-í) J 2k e dxj 



Sa Jq = ti-lfVïk* 



Wj 



ós 



1 H 



Jq 



R.J 



(1) 



with k^ làser wavenumber, wj spot size, Rj mirror cur- 
vature radius and 89 j q = (Sflj x x) ■ q depending on the 



rotation's angles SQj. In addition the mirror vibrations 
are accounted for by introducing matrices çj s and func- 
tions Sçjs representing respectively profile and amplitude 
of the s-th vibrational mode of the J-th mirror. 

The radiation pressure force and forques are linearized 
with respect to the set {Saj q } of transverse mirror 
coordinates by introducing optical stiffness coefficients. 
Hence, the Fourier transforms {Saj q } satisfy Langevin 
equations with driving forces proportional to these quan- 
tities. They become important in proximity of the sus- 
pensions and mirror mechanical resonances. These stiff- 
ness coefficients are in general frequency dependent for 
the presence of the phase factor e mT , with r the cavity 
round trip. 

A vector approach has been adopted by representing 
the amplitudes of the excited cavity modes by a column 
vector a while the mismatch and misalignment of the 
input beam is accounted for by a vector v. The vari- 
ous quantities O used for describing the system dynamics 
(e.g. the force acting on a given suspension mode) have 
been expressed in forms like O = v' ■ O • v with O a 
matrix representation of the quantity itself. In analogy 
with quantum mechanics O is seen as the matrix repre- 
sentation of an operator O corresponding to the quantity 
of interest, and v as the quantum state of the resonator. 

For evaluating spectral densities it has been introduced 
the symbol * which is defined by its action on frequency 
dependent quantities: 

/* M = /* (-w*) , 

and the shorthands 

9?{/H} = ^(/H-/ í M) 
K{/H} = ^(/H + /*M) . 

The same $ applied to a frequency dependent matrix 
transforms O (w) into O* (w) — (— zu*). 

The summation symbol is omitted when applied to ex- 
pressions containing a repeated index. 

The paper is organized in seven sections. Section lïïl is 
dedicated to the optical modes excited in a cavity with 
moving mirrors, and to the susceptivities relative to the 
noise sources of the suspensions, mirror vibrations and 
shot noise. The dynamic of the mechanical components 
(mirrors and suspensions) is discussed in Sec. IIIII while 
the Drever-Pound and quadrant detector signals are an- 
alyzed in Sec. IIVI 

The results obtained in these sections are combined in 
Sec.lVlwhere a five degrees of freedom model of the cav- 
ity, including radiative pressure and forques, is presented. 
The model is linearized for small misalignments and res- 
onance enhanced effects are discussed in Sec IVII where 
the cavity is examined as a bipartite system. In this 
context ponderomotive squeezing of the output field and 
entanglement of two mirror modes are discussed. The 
manuscript is completed by six mathematical appendices. 
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The first three of them give the expressions for the stiff- 
ness and the Drever-Pound signal matrices together with 
their simplified expressions in case of small misalignment 
and mismatch. The last ones are dedicated to thermal 
and shot-noise sources, and their mutual correlations. 



II. CAVITY FIELD 

A suspended cavity of length L excited by a time har- 
mònic field is described by a superposition of Hcrmitc- 
Gauss modes u\ (r, x) 



E 



exp 



iak 



2R(x) 



(x, t) u x (r, x) 



with oj e the làser frequency, a = + for a wave travel- 
ling from mirror 1 toward 2 and — contrari wise. The 
wavefront curvature R(x) is matched to the mirror's cur- 
vature: R(0) = ü x < 0, R(L) = R 2 > 0. Each 
mode is labeled as usual by a couple of integer num- 
bers (X y , X z ) = A. Here and in the following x stands for 
the optical axis coordinate and (y, z) = r are the two 
transverse ones. Each mode u\ is taken with a fixed nor- 
malization on the transverse section and without phase 
factors, 

u\{f) = u Xy {y)u Xz { z ) , 
while the amplitudes are written as 



» 



(x,t) 



^l<7\k e X-(\ y +\z + l)4>(x) 



» 



(x,t) 



where (x) = arctan ( x ~ b x ° ) is the phase delay of the 
Gaussian fundamental mode with respect to a plane 
wave, xq being the distance of the waist from the in- 
put mirror and b the confocal paramctcr. The field is 
propagated outside the resonator by passing through the 
different optical components met on the way toward the 
làser source and the photodetector which provides the 
error signal. 

The làser beam incident on (input) mirror 1 has been 
split it in a classical and in a quantum term: 

E m (r, i) = e'^E (l + // (*)) u m (r) + Sà SN (2) 
being 



E 



P Q f P 

^/= 2 · 5x10 (ïw 



1/2 



l f J 



1/2 



Hz 1 / 2 



the mean amplitude, and /i the relative amplitude fluc- 
tuations. Effects of the làser linewidths have been ig- 
nored. 

Misalignment and mismatch effects between the input 
beam and the cavity are taken into account by writing 
u ln (r) as a superposition of cavity modes, namely 



(f) = V\ u x (r) 



(3) 



The structure of the expansion coefficients is factorized 
in a product of Hcrmitc polynomials 



V\ oc 



S y S z H Xy (v„^) H Xz (v z ^ 



^y+^XylXJ 



(4) 



depending on the misalignment v q and mismatch 8 q pa- 
rameters defined respectively by 



fe W\ 



Q q 



5 q = Jl + i 



QqQl 



For a perfect matching v = and 6 = 0. Here Qi is the 
complex curvature radius of the cavity mode evaluated 
at the input mirror, while Q qi 9 q and e q stand for the 
curvature radius, angular and transverse misalignment 
of the input beam. 

The modal expansion J2Jl will be used in the following 
for representing the cavity fields in correspondence of the 
two mirrors as column vectors v with components v\ . So 
doing the multiplication of u (r) by a function w (f) will 
be represented by the product w • v of v by a matrix w. 

The coupling of the cavity with the universe modes 
through the partially transmitting mirrors introduces the 
quantum noise contribution 5à SN (r,í) of Eq. (J2J 

[Sa SN (r,í) , 6à SN ï (i-'/)] = <S (3) (r ~ r') S [t ~ i!) 

It can be expanded as a superposition of delta-correlated 
operators Sà\ N (t), 



6à SN (r,í)=íar (t)ux (r) 



SN 



(5) 



Before arriving at the mirror the excitation beam is 
passed through a phase modulator represented by the 
phase factor 



F = e lMsiniAt) = e- ipAt J p (Aí) 



(6) 



with J p (M) the p-th Bessel function of argument M. 
The input modulation F modifies the làser excited am- 
plitude a e x into a sum of harmonies varying on the time 
scale of the suspension fluctuations, 



(7) 



a\ = J P e lpAt a Xp 
while leaving the noise unaffected. 

A. Cavity fluctuations 



Owing to the fluctuations of the suspensions the mirror 
orientations change slowly in time by undergoing torsions 
Sfljz (í), tiltings SQjy (t) and transverse displacements 
5sj(t). The mirror can rotate also around the optical 
axis, but this motion is uncoupled to the cavity field in 
the linear approximation. 
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The mirror motions separate into fluctuating and av- 
erage components, the latter ones setting the reference 
frame for the vector representation. So doing the aver- 
age misalignment and displacements will be included in 
those relative to the input beam, which will be repre- 
sented by a unit amplitude vector vj 



vi 



with <fr a diagonal matrix of components = 
e - i 2(A B +A,+i)0 G and ^ (o) the single-trip 

phase delay of the Gaussian fundamental mode. Accord- 
ingly in the following the parameters 5ctj q (see Eq. Q 
will be small fluctuating quantities. 

The reflection at mirror f induces the transformation 
Uxdx = x: i u \ a \ r ' > w ith r i (t) the phase factor 



ti = n exp 



_ atfSxí - i2k l 5uf EF 

Ri 



m' ( süi x x - i / 



(8) 



r 



Here fei (í) is the deviation of the center from the po- 
sitions at rest {x\ (í) = + Sx\ (t)). Su^ EF (r, t) is the 
tiny deformation of the mirror surface represented by the 
matrix <5 ç x (i) of components 

Sç lxy (t) = 2fc £ J u x (r) Su EEF (r, t) u x , (r) d 2 f (9) 

Expanding further Su EEF into mirror modes [2^, 0] 
Sç\ (t) becomes a superposition 



Sçi (t) = Sçi s (í) ç ía 



(10) 



of matrices ç ls times fluctuating amplitudcs 5çi a (t) 
driven by radiation pressure and thermal noise. 

Although the freqüències of the mirror acoustic modes 
are very large, the tails of their spectra contribute to the 
low frequency thermal noise of the interferometers as re- 
cently reported by [2í|. Levin J26j has approximated, 
at very low frequency, the many mode profiles with the 
steady-state mirror surface deformation Su EEF (r) (in 
matrix form çf ) under the action of the incident beam 
(positive for a compression), by replacing Ea. ljTüj) with 



Sü (t) = çf íçf (í) 

Sç L (t) being a stochastic process, (see Ea.(|43fl) 



(11) 



Accordingly, ignoring the quadratic expression 
k l 5e\/Ri the phase factor ti (Eq.©) is represented in 
vector form by 



ne 



.Di(-dai) 



(12) 



with Sa± = (óaiy, Sai z ) the combination of rotation and 
displacemcnt dehned by Eq.QJ and Di the displaccmcnt 
operator 

Di (Sai) = exp (-ía lg Bj + Sa* lq B q ) 

acting on the functions of the transverse coordinates. 
The operators B y and B z act on the mode functions u\ 
as typical annihilation operators, Bit„ = y/ntí n _i, and 
this is the reason why D has been called a displacement 
operator. 

Next, the propagation from the input mirror to the 
opposite one is described by 



(13) 



with A — e~ T at the delay operator by the cavity round- 
trip time r. Next combining i|12|) with i|13[) a round trip 
is represented by 



Re 



-íljj—iSlf)!,, 



A* 



-i(5ç2 



D 2 (Sa 2 ) ( A* ) * • e~ iSçi • Ü! (-í ai 



r 



where tp is the detuning phase (tp > for a cavity shortcr 
than the closest resonance length), R — r\r-2 — e^^^ 
with T the cavity finesse, and 5ipi,cav the accumulated 
phase shift, positive for decreasing cavity length, 



Hl,cav (*) = Sípi (t - t) + Sifo t 



with Stpj 



(~1) J 2k i Sx J . Next, in view of the small- 



ness of Sifji^avjdaj and Sçj Di 2 and e~ lS '^ 1 · cav can be 
linearized thus obtaining for the round-trip transforma- 
tion 



e-**ira (pt - iX ■ 5aj, cav - iSçj,cav) . (14) 
Here X • Saj tCav indicates the sum Xi (óctj tCav ) i and two 
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vectors 



( 1 1 X a i X 2 , ~Y V , ~Y Z 



(15) 



5a JiCav = (S^ JtCavl Saj ycav ,5a[j zcavl Sa'j ycav ,5a' Jzcav ) 

collect the phase quadratures X g = B 9 + Bj, Y g = 
i (B ç — Bj) and the combinations 

<fai g , CCTO (í) = óaiq (t - t) + e l4 ' G 6a2 q (t - ~) 

(c/ and a" are the real and imaginary part of a respec- 
tively) . 

Analogously for 5çj, cav 



Sçi,cav (t) = ÇjsSçjs (í - T) + * 2 • Ç 2s ' • * 2 fe 



Here Saj is fluctuating with the cavity geometry and 
làser intensity, while ao,j does not depend on it and on 
shot noise, 



ao,j 



£G ■ v,F 



ta., =a 5 (//G-i<5 • Saj, cav -i5G DBF j -vjF 



6& 



SN 



tiG- íaj 



(20) 



Furthcr, the relation b t e~ lpM = e~ lpAt e ipAT D t im- 



plics Ge 



-jpAtQ^ su ffix p indicating 



that R has been replaced by R p = e lpAr R. Then, the 
factor e~ lpAt contained in the function F (see Eq. [üj can 
be displaced from the right to the left side of the above 
expressions by adding the suffix p to the various Green's 
functions. Hence 



The amplitude (t) of the A^A^-th mode is prop- 

agated back and forth the cavity. The fraction t\ is in- 
jected into the Fabry Perot through mirror 1 at time t, 
propagates toward and is reflected by mirror 2 at t + J 
and again by 1 at t. Hence, summing over the sequence 
of rund-trips, the held aj incident on the J-th mirror 
reads 



a./ = £ (1 + ,/) Gj-vjF + Ssl sn 
with £ = t\E and 

^ 1 



(16) 



1 - Re [Dt-iX- 6aj, cav - iSçj, ca v 



For very small Saj yCav and 6 Çj lCa „ first-order pertur- 
bation theory can be used. On the other hand assuming 
for u\ either Hermite or Laguerre-Gauss modes the var- 
ious terms of the perturbation X ■ 8aj yCa v + 5çj,cav do 
not couple the respective degenerate modes. Hence, the 
Green operator Gj can be expressed as 



Gj ~ G— i0 ■ Sajcav—iSG 



DEF 



(17) 



where the first term on the right is a static propagator, 
the second the contribution of the linearized motion of 
the mirrors and the third one describes the mirror defor- 
mat ions, 



(ao,j P + 8a.jp) + Sk 



-SN 



where 



a 0,Jp 



£J p G p -vj 



(21) 



(Sajp = e ipAt £J p (l/Gp — i<3 p ■ 5aj, cav —iSGp EF ^J -vj 
Analogously for the output field [Ï6| 



^OUT 



iOUT 



o,j p — t\£J p G p " -vi 

= h£J p (Gp UT f/ - Í0 P ■ 5a x , c 



í)a OUT 



-lòGp I -vi 



Sa 



OUT SN 



= t 



2 I qOUT.x-SN i t2 r<_xzSN 



(22) 



where G OUT = G - n/í?. 



III. RADIATION PRESSURE AND TORQUE 

Bouncing back and forth the two mirrors the làser and 
shot noise fields exert a radiation pressure resulting in an 
axial force directed along the òptic axis x and a torque 
parallel to their surfaces, proportional to the total inten- 
sity a^-aj and moments a^ • X ç • a,j. They split into 
classical Tj rp , Tj rp and quantum T^ p N , T§ £L compo- 
nents respectively given by 



SG 



G 

© 

DEF 



1 - Re-^Dt® 



e-^RG ■ <& X G 
e-^RG ■ $-óçj, cav ■ G 



(18) 



J rp 
Tj rp 



= (-1) J £ 2 2mj hk E (F ,j + SFj) x 
= (-l) J £ 2 2mjhk^(T 0Jq + ST 



and 



(23) 



Jq) qxx 



Next, the contributions of the shot noises entering the 
cavity through mirror J has been split as 5à SN = 5&f N + 
Í2Í 1 " 1 i5àf N , so that the same approximation of Eq. I|18|) 
applies and 



a./ = a ,j + Saj + 5a 



SN 



(19) 



<pSN 
J rp 

n-SN 
J rp 



{-l) J £2ÜKj Hk e X^x 
(-l) J £2mjhk^X s je N q qxx 



(24) 



where 9tj = | z',/ 1 2 + \ Aj with Aj the J-th mirror power 
absorption. 
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and aj , 



i<o,j and SFj indicate the contributions of aj j-slq.j 
S&j + H.c. and analogously for To.j q , 5Tj q . 
Fo,j,To^j q split in turn into time constant terms Fq, 
To,,/ ? , balanced by the stabilization system of the appara- 
tus, and small terms SFq^j, STo,j q oscillating at múltiples 
of A. Being A typically of the order of some MHz these 
contributions can be ignored. 

For a stabilized resonator Gj is represented as in itTTI) 
so that SFj and STj q reduce in the frequency domain 
respectively to 



SFj = F ,jtf 

ÒT Jq — T 0t j q fl 



$j ■ Sa J<cav + 5F Jcav 



DEF 



•-.i,, oàj. cav + STj q E c F v 



(25) 



The three pieces of Eqs. (I25|) represent, in the given or- 
der, the contributions of the fluctuation of làser intensity, 
mirror displacements, rotations and surface deformations 
to the radiation pressure forces and torques. 

Being the suspension characteristic freqüències gener- 
ally smaller than the mirror modes resonances [24| . the 
deformations are described by a single matrix (see Eq. 

GD)- 

The vectors $j = (Fj^, F JXq , F JYq ) and Í Jq = 
{Tjip q , Tj q x q > , Tj q Y q > ) contain five proportionality con- 
stants between the forces (the torques) and the coordi- 
nates (Sàj tCav ) which parametrize the mirror's displace- 
ment, so they are stiffness coefficients. F _j, SFf^ 

and To,jq, Íj q , 5Tj q E c F v depend on the steady-state am- 
plitudes of the cavity modes, represented by the vector 



O 0J 
Òj 



= vt 



O • vj 
Ò vj 



(26) 



SF? c E a F , T Jq and 



with O 0: ,/ = F 0:J , T j q and O = $j 
ST? q E f av Matrices F ,SF EEE and T 5T% EF V _ are re- 
ported in Appendix A (Eqs. (|A1IA4|I ) while Í q are 
collections of five matrices (Eq. (|A2|0 . They depend on 
Green's matrices (Eqs. (|A3IA5|) 'l. and through them on 
frequency and detuning, closeness of cavity modes with 
respect to linewidth and phase modulation dcpth. The 
frequency dependence is due to the factor e mT appearing 
in different fashions in G p , Gp UT , G p . 

Eventually, the shot-noise contributions (Eq. (|23}) are 
expressed by 



Xj^ — í 1 1 JpVj • G^-Xj • S 



gSN e i P M 



H.c. 



with i € (ip,6 y ,6 z ) and take in the frequency domain the 
form 



vSJV _ ±- 
^Ji — r i 



tr 1 2J í ,3?{v^·Gt·X i ·^ JV } 
Finally, on the J's mirror mode act the forces, 



(27) 



T° EE = (-1) J 2Kj hk* 



F 



DEF 
0,Js 



SF 



-pDEF SN 
*' Js rp 



(-1) J £2dK J hk í Xf EF SN 



DEF \ 
Js,cav } 

X 



(28) 





L(cm) 





FIG. 2: Axial F^,, F 1Xy , Fixz Fxy v and angular T lzXy , 
TizXz, TizYy, TizYz stiffness coefficients vs. length L of a 
symmetric cavity for angular misalignments 9 y — .01, #2 = .1 
mrad and detunings tp = .1, .2, .3 -n . The round-trip phase 
factor e ïtCT has been ignored. 



where 



xpDEF 

u Js,cav 



+ é 



r/2 fpDEF c~_ 



Hcrc 



ÒDEF 
JsJ's' 



v t .fDL·l 



JsJ's' ' V J 



(29) 



is the force acting on the Js-mode due to the deforma- 
tions of the mirror surfaces. In this case the force does 
not factorize as for the suspension modes. FjfjT s i (Eq. 
(IA6I) ) represent the effects of the vibrations of the modes 
J's' on the Js one. 

Next, the shot-noise force is given by 



yDEF sn 
A Js 



t^2J p lSt\w 



GUjs 



<5à. 



~SN 



} (30) 



In Fig. |21 the optically induced stiffness coefficients 
have been plotted for a set of detunings and angular mis- 
alignments, in an almost concèntric cavity having a fi- 



nesse T = 500 and output spot sizes of 2 x 10 



Be- 



ing close to the concèntric configuration also the stiffness 
\xy,T y Yz become comparable with 

10~ 2 rad, 



coefficients F x /Y q jT i 



FipiTj q x/Y q for cavity axis misaligned by 6 y 
6, = 10 _1 rad. The signs of the stiffness coefficients may 
have important consequences on the mechanical stabil- 
ity, as discussed by several authors for plane-parallel and 
concave mirrors [ÏCl Í18| . 



A. Small misalignment and mismatch 

In thc limit of small misalignment and mis- 
match the vector Vi (Eq. (0J) reduces to ~ 
1+Sv ={1, (Sv y ,8v z ),0, ...,0} with 



S(a.u.) 



V2 Q*Q q 
Wl Q* - Qq 



Ri 



For mirror 2 v ç is multiplied by e _í *' G . Splitting forces 
and torques in 0-th and 1-st order terms in these mis- 
alignment parameters SFj and 8Tj q of Eqs. (|25|l take 
the simpler forms, 



5Fj 

sf Jq 



5T 



(o) 



Jq 



ST 



d) 

Jq 



(31) 



where 



SFj 0) = F ,jF + Fi l 61>j, caa +5F%£Z 
SF ( j 1} = 5F JXq Sà% 

,cav + 5FjYq5à' Jq 

,cav 

i rp c ~ I 

l Jq,cuv ' Y^^Jq,cav 



DEF 



6T™ = fxsm 



with 



SFjx/vq = 2Re|F x:/ yVj g | 
Sfj q = 2 Re [t^vjgj 

Fx/YiT<p being defined in Appendix C. Accordingly, in 
the ideal setting of the cavity the forces and torques 
are respectively proportional to longitudinal 8ipj Mav and 
transverse Sà'j g cav fluctuations through the stiffness co- 
eiïicients F^,, T x /y- A slight deviation from it introduces 
forces and torques with a reverse dependence on fluctu- 
ations, say SFj 1 ^ , 6Tj£ depend respectively on 5àj q . cav 
and Sip^ cav . 

IV. ERROR SIGNALS 

The errors used for controlling the cavity are provided 
by Drcver-Pound (DP) and quadrant detector signals 
(QD). In thc DP detection technique the photodetector 
current I(t), obtained from the light transmitted and 
reflectcd by the input mirror, is mixed with a local os- 
cillator ~ sin (kAt + tp) with positive odd integer k and 
low-pass filtered by an averaging procedure 



-.DP 



(*) 



K DF (t - t') sin (feAí' + ip)I (t') dt' (32) 



with the filter response K DP (t — t') extended to a suit- 
able interval much longer than (fcA) , and short com- 
pared to the time scale of the phase-quadrature fluctua- 
tions. Tuning (p around s DP can be maximized for a 
misaligned cavity. 




FIG. 3: Drever-Pound static characteristics s vs. tp for cav- 
ity lengths 1 4- 10 cm and (p = 0. The plots correspond to 
modulation frequency A = 2c\ e / [nw 2 ), depth M = 0.1, and 
jaw angle S9 Z — 0.01. 



Putting 5&p — 5íi b (w +pA), s is represented 



the frequency domain by 



~DP i Í>DP 



where 



S5 DP DEF 
j^DP SN 



= vt-I 

= vï 



DP 



zDP 



ss 

+£X DPSN 



~DP DEF\ 



(33) 



• vi 

ínDP ^DP 



si DP DEF 



X 



DP SN 
+ 



- j: ' ) • vi 

_ fifDP DEF 

DP SN 



Vl 



X 



with I DP ,SÏ^ PDEF ,3 DP ,Xg PSN defined in Ap- 
pendix B. 

In Figure[3]the static characteristic s DP versus ip has 
been plotted for a set of cavity lengths and modulations. 

Figure 4 contains plots of the coefficients s$ and sxy 
vs. cavity length for (p — 0, 9z = .1 mrad and 7 de- 
tunings. They show that as a consequence of the mis- 
alignment sxy becomes comparable to s^, so that the 
D-P error signal contains contributions of the torsional 
fluctuations around the vertical axis. 

At low frequency Ss DP DEF becomes proportional to 
the thermal noises (5ç L TH . 

The quadrant detector used for stabilizing thc angular 
oscillations provides two error sig nals sg (í) (q = y,z), 
proportional to expressions similar to (|32|l with the cur- 
rent / replaced by I q — /3' ■ Q q ■ [3 with 

(i = e-^ At a^ UT + 6k OUT SN , 

the matrix Q 9 representing the function sgn (q) . Thcn, 
D is given by an expression similar to (|33|l with Gp UT ' i 
replaced by G° C/Tt Q q . 



A. Small misalignment and mismatch 

In the limit of small misalignment and mismatch the 
signal can be split into zeroth- and first-order contribu- 
tions 



-DP 
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FIG. 4: Coefïicients and sx y vs. cavity length for ip = 0, 
#2 = .1 mrad and 7 detunings ranging in the interval —.3-^ < 
V<-3§. 

given respectively by (Eqs. I|()4|l) 

g(0)/^P = £ 2 + ~ s ^ l<cav + S ~ s (0) DEF^j 

+£X SN W 

SsU/K DP = £ 2 (6s Xq 5à'^ cav + 6s Yg ó&[ qtCav 

+ S§W DEF } + £SX SN W (34) 

with <5sx/y q = 2Re {sjc/y^ig} an d defined in Ap- 

pendix C. 

For a perfectly aligned and matched cavity the D-P sig- 
nal is sensitive to the axial fluctuations 5ipi lCav , mirror 
deformation term Ss^ ' DEF and shot noise X SN (°> . In 
particular, a contribution Sipx tCav depending on the mir- 
ror thermal noise Sç[ 2 is added to the length fluctuation. 
A deviation from alignment introduces in the error signal 
contributions proportional to the transverse fluctuations. 

V. 3D MODEL 

The deviations of each mirror from the reference posi- 
tion is described by the displacements Sx, Se y and 5e z of 
its vèrtex and the angular parameters S9 Z = —5íl y and 
66 y = Sft z . As said Sfl q describes a right-handed rota- 
tions around the axis " q" , so that S9 Z is a left-handed tilt 
and 59 y a right-handed torsion. These quantities fluctu- 
ate as a consequence of suspension thermal fluctuations 
and mirror surface deformations. The radiation pressure 
fluctuations are transferred to the mirrors proportionally 
with the làser intensity. The cavity reacts by changes 
of geometry which in turn changes the stored field and 
closes the loop of the cavity-field system. 



From a purely-mechanical point of view if the design is 
good (that is, symmetric enough) the suspension masses 
are aligned along the vertical axis z, perpendicular to 
the cavity axis x. In thcsc conditions the torsion S9 y 
and vertical 8e z degrees of freedom are uncoupled. A 
coupling between longitudinal motion Sx and tilt S9 Z is 
generally speaking unavoidable. This is true also for the 
transverse displacement Se y which is coupled with the 
rotation around the optical axis. It goes without saying 
that in a real situation it is very difflcult to avoid more 
general cross couplings. 

Radiation pressure can increase or reduce existing cou- 
plings, and it can also produce new ones. While Se q is 
insensitive to radiation pressure, S9 q responds to the ra- 
diation torque. For this reason when asymmetric optical 
modes are excited the rotations 59 q modify the radiation 
pressure, and ultimately couple Sx and tilting, but also 
Sx and torsion. 

Before proceeding further it is worth replacing the dis- 
placements 5ej q by 6ipj q = 2k Sej q , the angles S9j q 
by S'djq = \/2k í wj59j q and introducing a new five 
component vector Sípj = (Sipj,Sipj q ,S'dj q ) which forms 
with the cavity mode amplitudes a system of correlated 
stochastic processes. It is usually a very good approxima- 
tion to model the suspension as a set of damped, indepen- 
dent oscillators coupled to an heat bath. Each oscillator 
JA j, labelled by A, specifying the prevalent character of 
the mode (tilting, torsion, displacements, violin modes), 
and the mode index j, can be parameterized with its effec- 
tive mass Mj\j, pulsation vjj\j and damping coefficicnt 
lJ\y For rotations Mj\j is replaced by the moment of 
inertia. These parameters are related to the masses and 
stiffness constants of the system. The coordinates of the 
mirror can be written as linear combinations of the oscil- 
lator 's coordinates q n , and this means that each normal 
mode gives in principle a contribution to the mirror's mo- 
tion. By interacting with thermal baths these modes un- 
dergo Brownian motions by influencing the electromag- 
nètic field, eventually coupling mechanical and radiation 
pressure fluctuations. 

A. Suspension Langevin system 

By linearizing the equation of motion of each mirror 
(J) the horizontal (x and y) and vertical (z) displace- 
ments Sípj, torsion Siïjy, tilt S~dj z and rotation around 
the cavity axis ó$j x are expressed in terms of the ampli- 
tudes Ajxj of the normal modes as 

Sípjp. = K Jl·l \jAjxj 

having indicated by J/i a genèric degree of freedom and 
by Kj^xj the coupling coefflcient with the mode JAj [23] . 

If the mirror vèrtex coincides with the center of mass 
of the suspension payload, and the centers of the sus- 
pended masses are aligned along the vertical z-axis, the 
suspension can be easily modeled by considering only the 
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couplings 5ipj — S'Sjz and 5ipj y — S'djx, and assuming the 
vertical oscillations independent of the other degrees of 
freedom. Being the amplitudes of the cavity modes in- 
dependent of the rotations Sií j x , the suspended cavity is 
described by the collection Sípj of five fiuctuating quan- 
tities, depending linearly on radiation pressure-torques, 
thermal noise, D-P and quadrant detector error signals, 



Hj = XJï,8mj(£ 2 6Fj+£Xj j 

?2<;rf, i CYSN\ 



+Xj^ + S J1 H DP S DP 



TH 



Hjy/z = X Jy/z 

5$j z = Xje zi) 8mj(s 2 SFj+£X 



SN 



+XJ9?mj [£ z 8T Jz + £X 



JBz J 



60 



Jy 



(£ 2 Sf Jy + £X 



SN 
JOy 



(35) 



In writing Eq. (|35(l the interaction with the mirror 
noise was approximated with Eq. while in Eq. I|3t)|) 
the effects of the suspension fluctuations were ignored. 
Loosely speaking the two systems refer respectively to 
the low and high frequency regions. In the former one 
the suspensions are mutually coupled by radiative forces 
represented while the mirror vibrations generate a global 
thermal noise hiding the single mode contributions. In 
the latter one the suspensions appear frozen and the mir- 
ror modes are mutually coupled by radiative forces rep- 
resented by 6F? a ?£ v . ' 

The solutions of the homogeneous system <|35ll repre- 
scnt, in absence of feedback forces, free mechanical os- 
cillations of the suspended cavity, stable or unstable in 
accordance with the sign of the imaginary part of the 
oscillation frequency [2(j ■ 

For a more detailed analysis l|35|l and (|36[1 should be 
mirrored by the system relative to the quantities í^j > 
óçj conjugate of 5tpj, {6çj s } , which can be obtained 
from the above one by replacing x.Ji>/e q / s by X%/e q / a 

(Eq. m and XTB^xDBPTH by yTH^yUEFTH 

(Eq. I41fl in the random force expressions. 



In case the mirror vèrtex and/or the centers of the sus- 
pension wire clampings are displaced from the respective 
mass centers, the vertical fluctuations are coupled to the 
other ones. 

The effect of the servo systems acting on the longitudi- 
nal and angular mirror displacements have been included 
by indicating by H DP and H^ D the respective transfer 
functions. 

For the mirror vibrations a Langevin equation for each 
mode must be considered since their profiles are different 
(Eq. ESJ), 



6? Js = 8xjs (£ 2 5Ff s % u +£X 



DEF SN 
Js 



vDEF TH 
<*Js 



(36) 

Expressing 5Fj, 8Tj q in terms of displacements and ro- 
tations by introducing the stiffness coefficients, and doing 
the same for the error signals s DP , s® D the above system 
can be reduced to an equivalent one relating the fiuctu- 
ating displacement + rotations to the thermal noise and 
shot noise sources. 

The axial displacement Sípj and tilting ó·djz re- 
spond to the axial force £ 2 SFj + £Xj N and torque 

( £ 2 5Tj z + £Xj^j y generated by the làser beam and 

the shot noise. By the way they include the contribu- 
tions of the mirror thermal noise. On the other hand, 
Sèjy responds to the torque £ 2 STj y + £XjqL. The links 

between force-torques and Sípj are represented by the 
susceptibilities xj m „. 

The terms proportional to £ 2 and £ describe the re- 
sponse of the system to radiation pressure. Their pres- 
ence indicates that a motion of the mirrors produces not 
only a phase change but also an intensity change provid- 
ing a spring action. 



B. Susceptibilities 

The susceptibility xj^ v describes the action on the co- 
ordinate /i of the force/torque acting on 



XJfiv 



with x,j\j the susceptibility of the mode JXj of frequency 
and damping coeficient 7,7 aj 



XJ\j 



_ „LD 2 



(37) 



and Kj^xj, Kj„xj the coupling coefficients with \i and v 
mirror coordinates, while the adimensional Lamb-Dicke 
factor 



LD_ k i 



VJX3 



ZMjxjZVjxj 



(38) 



depends on the mode mass Mjxj = MjiK 2 JiX ~ (the subfix 
i idcntffies the i— th mass of the suspension) . For rotations 
Mjxj is replaccd by Ijxj/w 2 with Ijxj the moment of 
inertia. Some authors use the so-callcd optomechanical 
coupling constants Gjxj — 2v / 2?7jfj/" r pïj . 

The mechanical susceptibility x.Js is similar to 1)3 7|) 
while the mass appearing in the Lamb-Dickc factor varies 
for the different modes, as reported in j24[. 



C. Thermal contributions 

Assuming suspension masses at the same temperature 
T, each mode is characterized by a thermal source (see 
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Appendix D) 



vTH _ 



Ak B T 



JXj 



Oxj 



IJjXj 



?JAj 



3k B T 



JXj - 
VJXj 



(39) 



with fjjí delta correlated random forces introduced by 
Diosi |12j in order to remove some inconsistències of the 
classical Langevin equation. 

A Y-version of ll-iól) can be easily obtained for the Y- 
quadratures corresponding to the above ones by replacing 
XJ^Xj by 



XjXj 



_ 

l X.JXj 



(40) 



For a Gaussian illumination cp takes the form 

\-o 2 

CG : 



2nEw,j 

with wj the spot-size and E,a the Young's modulus and 
Poisson ratio respectively. For a genèric illumination Cp 
can be expressed as cp — fp cq with 



, P: 



A„,A Pa' A 



XX' 



p2 



A„ + A' A,+A' 



(45) 



/q/3 being the af3 cocfficicnt of the expansion of 
5ug (r) e - ^- m moc i es u\ \ t (r). 



and X T j* by 



vTH 
ï JXj 



1 4k B T - 



'JAj 



3fc^T í?JAi 



(41) 



The terms proportional to ryjj in Eqs. i|39|) and 1411) 

can be generally neglected except when the temperature 
is rather low and the oscillation freqüències very high, a 
situation met only in some mirror modes. 

ijjxj disappears in the simple Brownian motion model 

while in Ref. Q fjjxj has been dropped and J ~~^j^À JXj 
replaced by a new delta correlated random noise source 
Qjxj- 

The thermal sources Xj^ are superpositions 



X 



TH 



L J M ~ -^J/íAj XjXj X JXj 

of the Xjjl weighted by the thermal susceptivities 



TH yTH 

v.-A 



XjXj = K JAj XJAi 



(42) 



With Kjxj = Zy/ÏJXjhjXy 

The terms of (|35|) contain contributions proportional 
to the fluctuating quantities 8çj ™ |2(| 



4k B T 



\/2hk e cp(/)jç 



(43) 



with (f> the loss angle, ç a delta correlated random force 
and cp depending on the illumination profile 



P (f) = P; 



A„,A,e 



! "A„,A 



(r) 



(44) 



For P (r) differing notably from the Gaussian one the 
deformed profile of the mirror óuf.f F can be expressed, 
neglecting the finite size of the mirrors, by a suitable 
combination of derivatives of the deformation ÓUq EF (f) 
relative to a Gaussian distribution 



8u 



DEF , 



£ 

Aa.Aj. 



X y ,X Z 



-w) 



X y + X z 



dy x y dz Xz 



VI. THE SUSPENDED CAVITY AS A 
BIPARTITE SYSTEM 

When the frequeney is in proximity of two close reso- 
nances of the mirror 1 and 2 modes, the system behaves 
as a quantum mechanieal bipartite system described by 
Gaussian continuous variables. These systems can form 
EPR states characterized by their covariance matrix a 
which can be used for evaluating the entanglment of the 
state and its content of quantum information. 

The difference between the e.m. fields used in quan- 
tum òptics and the present mechanieal system concerns 
the sources of the respective states. The e.m. fields are 
produced by the e.m. vacuum noise entering through 
the mirrors of a cavity containing a nonlinear crystal. In 
the present case thermal and shot noises act as sources. 
Accordingly, the covariance matrix a can be split into 
thermal a TH and shot noise (8£) 2 a SN contributions ob- 
tained by separating Sçj into 5çj — 8£Sçj N + 5çJ H sat- 
isfying the Langevin system (j3*üjl 



- 5 ~SN/TH- 


1 


g -SN/TH 


~~ D 



V22 
-V21 



-V12 

V11 



XiX 

X2X. 



DEF SN/TH 
1 

DEF SN/TH 



(46) 

with Vjj> factors representing the radiation pressure ef- 
fects 

Pjj = 1 - 8e mT £ 2 Kj xjF?j EF 
Vjj = 8e m ^ 2 £ 2 m JX jF?j EF 



and their product D = VxyPii — VviVn- An analogous 
system holds for SçJ with \ J replaced by xS ■ 

' ~T23) 



The output field contains a component (Eqs. 



-OUT 



proportional to <5çi.2 through the matrices Zj = 
G ■ $ ■ ÇJ ■ G and a shot noise G OUT -6àf N ->- ííri-X&SN 



ti z 



term. Hence, depending <5çf2 linearly on the quadra- 



tures Xfi 



DEF SN 



the output exhibits some degree of 
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squeezing., a feature exploited by several groups in the 
context of gravitational antennas of the next generation 
[2^ | . The dependence of the efficiency of the ponderomo- 
tive squeezing on the mirror deformation profiles (ma- 
trices Zj) and residual misalignment/mismatch can be 
easily analyzed by means of Eqs. (|46|l and the correla- 
tions of Apps. E and F. 

The complex dynamics of cavity field and ponderomo- 
tive effects may lead to the creation of quantum entangled 
states of the two mirror modes, as shown by Mancini et 
al. ([IJ and references therein included). These authors 
have proposed a measure E (w) of the entanglement de- 
gree (the smaller E (zu) < 1 the larger the entanglement) 
based on a combination of the elements of the covariance 
matrix, 



EO) 



~Y 



(47) 



Splitting the quadratures into shot noise and thermal 
contributions, taking into account the many modes of the 
cavity and the shapes of the mirror mechanical modes, 
and scaling the ratio terms by keeping constant E(nj), 
yield for the thermal and the shot noise contributions 



\8^ H 


+ 6£f H \ 2 











= Xj 



TH | 2 (jTH X(-\ 



J 



(48) 



— J \a.jX T j HÍ 

SN X(+) 



[6ç? N ,SçÍ 



~SNY] 



XJXj'Cjji 

zo(zvj + mjt) _ _„ *^SNY 
-XJ Xj> a J <*i c h' 



Ivjjwji 



~Y TH 



and 



with Xj H defined in while|(5ç^ 
\SçY SN — dÇrf SN \ 2 are similar to (|48() -a and -c with xj, 

Ò\ 
C 



TH X(+) 
J 

TH Y(-) 
J 



and C 



SN X(+) 



and 



jj, replaced respectively by \j 1 

SN X{-) 



•n(±) 



(ai,a 2 ) = \Vi2,V\2 
Vjj ± Vjj and 



On the other hand, 

>(+)-5(+)<5(-)-5(-) _1/2 



C 



TH X/Y(±) 
J 



Re{óf X/YYTH } 



^SN X(±) 



r<SN Y 



Re 



;(±)*^(±) 



Im 



V), 



2J P^P^V^V^ 



with Cj^j, s , given by Eq. I|F2(1 . In App. E thermal noise 
correlations for the Lindblad-Diosi and the Giovanetti- 
Vitali MEs are explicitly given. 



VII. CONCLUSIONS 

A suspended cavity illuminated by a làser beam has 
been described as the mechanical response Sípj of each 
mirror of a linear system to radiative, thermal and shot 
noise forces. These perturbations have been linked to 
the mechanical responses by means of susceptibility co- 
efficients. The model includes the mirror vibrations de- 
scribed by a set of mode amplitudes Sç, together with 
their shapes ç. 

The radiative pressure forces and torques have been 
linearized with respect to Sípj and <5ç, by obtaining sets 
of stiffness coefhcients for the suspension and 1) and 
for the mirror modes {Ff^f a , ) . Accordingly the radiative 
forces have been expressed as produets of susceptibility 
cocfficients, làser intensity transmitted to the cavity (£ 2 ), 
stiffness coefhcients, and mechanical mode amplitudes. 
So doing Sípj and 6ç have been linked directly to the 
thermal contributions, modelled by Lindblad and non- 
Lindblad master equations, and to the shot noise forces. 
The mirror thermal noise has been expressed in the low 
frequeney limit by the Levin's formula. A corrective fac- 
tor, for taking into account deviations of the cavity field 
from the fundamental mode, has been introduced. 

The Drever-Pound and quadrant detector signals used 
for stabilizing respectively longitudinally and angularly 
the cavity, have been expressed in a form suitable to 
study the mutual coupling of these degrees of freedom 
in case of misalignment. 

Emphasis has been put on the description of missalign- 
ment and mismatch of the input lascr bcam. To this end 
a vector approach has been adopted: the state of the in- 
put beam and the amplitudes of the excited cavity modes 
have been represented by vectors (v and a, respectively) 
and all the contributions to the cavity dynamic by a set 
of matrices. In this way, all the relevant quantities are 
given in form of algebraic produets. 

In particular, the optically-induced stiffness cocffi- 
cients relative to the suspension modes have been ex- 
pressed in the form v^-Ç • v, v^-T • v with T matrices. 
It has been shown numerically that these coefficients may 
become very large in misaligned cavities close to unstable 
configurations. 

The finite cavity round-trip time has been included in 
the model by introducing a delay operator. Conscquently 
the cavity has been represented in the frequeney domain 
by frequeney dependent matrices containing stiffness co- 
efficients. 

The reported model simplifies notably in proximity of 
mechanical resonances. In particular the covariance ma- 
trix a of two close in frequeney vibrational modes has 
been expressed in terms of the stiffness coefficients and 
used for evaluating the system entanglement .This matrix 
also controls the squeezing degree of the output field. 

The numerical examples refer to almost concèntric cav- 
ities of length varying between 1 cm and 10 cm, spot 
size.2 cm and misalignment of .1 mrad. 
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APPENDIX A: FORCE, TORQUES AND 
STIFFNESS OPERATORS 



APPENDIX C: SMALL MISALIGNMENT AND 
MISMATCH 



F = JpGj, • G p 
To,g = JpGp • Xg • G p 

Ncxt, the stiffness operators % q are given by 



Assuming i? = 1 and p = the ponderomotive forcc 
and torques read 

|2 



(Al) 



F 



g, 



0(00,00) | 



2J, 



p 3 |t 9P | 



(A2) 



To,j q — 2Re |Go (00 00 -,G ( 10!l0 )U,/q| 
while the stiffness vectors reduce to 

Í, Jq - T(°)+2Re{^ 1) } 



(Cl) 



with 



while 



with 



$ p = e-^RpGl ■ G p ■ * • X ■ G p 
Í qp = e-^RpG\ ■ X g • G p • * • X ■ G p (A3) 



<5F 



DEF 

J : cav 

DEF 



2J 2 3 



Í?t?DEF \ 
I Jp.cav j 



ÍTjftl = 2J p 2 3{5f^} (A4) 



3 (0) = (^,0,0,0,0) 

= (o,FxVjy,FxVj z ,FYVjy,F Y vj z ^ 
= (f 4 {v Jv +v Jx ), 0,0,0,6) 



(C2) 



where 



with 



à^jp/av — e l ^RpG P G p <& <5ç^ cat , • Gp 
^T^ F caí) = e-^iZpGt X, G p * 5# • G p (A5) 

where Sçj cav takes the Levin's form, 

Finally the action of the modes J's' on the Js one is 
represented by the ensemblc of matrices 



2^ je ^ í2 '* G G' ( 00:00 ) | Co(oo,oo) | } 
23{e- i,/ '^ G G 0(10 , 10) |Go ( oo,oo)| 2 } 



Similar expression holds for Fx and with |Go(oo,oo)| 
replaced by G^ 00 00 \Go(io,io) while Fy is similar to Fx 
with 3 replaced by 9ï. 

Analogously, for the Drever-Pound error signal 



~ 5 DP 



s -(0) 

s (0) + 2Re 



{*>} 



FS5Í.-2JJ9{fJSJ , v } (A6) Js -dpz.ef = Js -mDi!F +2Re | s mi>EFj (C3) 



with 



V?f.j>.> = e" # ^Gt • ç Js • Gp • * • çj, s , ■ Gp (A7) 



APPENDIX B: DREVER-POUND SIGNAL 



with 

g(0) 
g(0) 

g(D 

íg (0) DEF 
SS (1) DEF 



(5^,0,0,0,0) 
(5^,0,0,0,0) 

(0, Sx^ly, Sx^lz, Sy^iy, SyWi z ) 

s 4> [ e Çis(oo,oo)0íis + e 



(C4) 



Ç2s(00,00)^2s 



1 DP = 2J p _ k J p Z {e^G p OÜT t . G out } 



{e-WTrt^GOt/rt . G p±fe * X Gp ±fe } 



SI 



DP DEF 



— *2<Jp±kJp 



p SÍV - 2Jp5R {ví • Gp UT ^ ■ SaÇïf SN ] (Bl) 



S X ( e ^ T ( w lyÇls(00,01) + S X Vl z Çls(00,10)) Sçu 
+é UT/2 («1^(00,01) + S X Wl z Ç 2s (00,10)) fes) 

where 

5 — f í oc> f „i·pr< OUT * r< OUT X 
*V — •Jp+k·Jp^ yi Lj p+fe(oo,oo) LT p(oo,oo) j 

S-0 — 2 t/p+ & Jptfi. |e ^ ^ G 



(e^iïp +fc G 



OUT * f'OUT r OUT 
p(00,00) ^+£(00,00) ^+£(00,00) 



' ï? C 
ri rl L·r 



OUT 



f" 



OUT QOUT 



P^p+feíOO.OO) ^p(00,00) ^p(00,00) 



)}(C5) 
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s x is similar to with G°^ 0000) , G^J m) replaced by 

G p+^io,io)' G p(ïcuo)' while S Y is similar t0 sx with 5i 
replaced by 3. 

Finally the shot noise contribution reads 

jflP SÍV _ jçDP SN (0) , j^DP SN (1) 



while for the master equation of Ref. [14j 
fixx th _ r< YY TH — r< XY TH 

On the other hand the commutators coincide 



where 



vDP SN (0) _ OT1fíír OUT * ÍÀn SN à?, sn \\ 

a — zj p jí |L·r p( · 00!00) ^oa p+fc ( 00 ) - oa p _ fe ( 00 ) j j- 

yflPSÍV(l) _ 9SÍ? Í^OC/T* ííx^SN ?~SN 
A — ZJt |«J p (io,io) \^ 0a p+fc(10) — oa p-k(10) 

+ (sà S p ? K01) -5~a s p » km )v* lz )} 



XDEF TH X /Y DEF TH 



= 4— (E2) 



APPENDIX F: SHOT NOISE CORRELATIONS 

The Fourier transforms of the shot noise force and 
torque Ij27|) are characterized by the correlations 



APPENDIX D: THERMAL AND SHOT-NOISE 
SOURCES 



1 + |)^ V <5(-W) 



X 



TH/SN 



V TH/SN - v 

+5 1J H DP £ 2 S DPTH / SN 



Jtl> 



TH/SN 



yTH _ - yTH 



X 



TH/SN 



yTH/ SN . ~ v ±n/ 
J8z - XJfl*V A J^ +XJBz^ J0z 

+5 lJ Hf D £ 2 ^ DTH / SN 
Xje H y = XJ*yXÏÏ + SuH? D í?$ DTa (Dl) 

while for the mirror modes 

S,DEF TH/SN V DEF TH/SN 

*Js = XJs^j s 



r TH/SN 



APPENDIX E: THERMAL NOISE 
CORRELATIONS 



The correlations X™ (w) X™ [w') 
c xx th s + w i^ _ of the thermal sources (O 

and (|4*T)l for the Diosi master equation (see 01) are 
given by 



çXXTH _ 4k B T t \zu + i-/j\ Hzuj 



W 



zu 2 r 3k B T 



(El) 



+2 



C 



YY 



th _ ^k B T ^ w 2 hzuj ^ ^txjj 
hzuj vj 2 3k B T vj 

M 



(jXY TH _ çYXTH*/_\ ^k B T 



tlVJj 



w 



3k B T 



where 



with 



f SN _ T , c sn . 



fiSN _ j2(~<* v 



G, 



and i, i' = 1,2,3. In particular, Cjfj ri , = Cj^/j^. 



Analogously for X 



DEF SN 
Js 



Xff^N {w)X DEFSN {wl) 



where 



with 



f<DEF SN 



+ fiDEF SN 



fiDEF SN 

WsJ's' 



J p ^p ' ^Js 

•Çj's' ' Gp 



G, 



In addition, 



XDEF SN ^) jX y,f, F SN {m>) 



í fiDEF SN\ r / , i\ 

\ c JsJ's> Jà{m + zu) 



(FI) 



(F2) 
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